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Any symmetry reduces a second-order differential equation to a first-order equation: variational 
symmetries of the action (exemplified by central field dynamics) lead to conservation laws, but 
symmetries of only the equations of motion (exemplified by scale- invariant hydrostatics), yield first- 
order non-conservation laws between invariants. We obtain these conservation laws by extending 
Noether's Theorem to non- variational symmetries, and present a variational formulation of spherical 
adiabatic hydrostatics. For scale-invariant hydrostatics, we directly recover all the published prop- 
erties of polytropes and define a core radius, a new measure of mass concentration in polytropes of 
index n. The Emden solutions (regular solutions of the Lane-Emden equation) are finally obtained, 
along with useful approximations. An appendix discusses the special n — 3 polytrope, emphasizing 
how the same mechanical structure allows different thermostatic structures in relativistic degenerate 
white dwarfs and and zero age main sequence stars. 

PACS numbers: 45.20.Jj, 45.50.-j, 47.10.A-, 47.10.ab, 47.10.Df, 95.30.Lz, 97.10.Cv 
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I. SYMMETRY REDUCES THE ORDER OF A DIFFERENTIAL EQUATION 

Noether's Theorem relates every variational symmetry to a conservation law, a first integral of the equations of 
motion, which can then be integrated directly by quadrature By an extension of Noether's Theorem, non- 
variational symmetries of the equations of motion also reduce them to reduced equations, which are not conservation 
laws Consider any system described by the Lagrangian L{t,qi,qi) and action S — J C{t,qi,qi)dt, where the qi 
are the coordinates, the dot designates the partial derivative d/dt with respect to the independent variable, and 
the Einstein summation convention is assumed for repeated indices. Under any infinitesimal point transformation 
6{t, qi), Sqj{t, qi) generated by dt ■ d/dt + dqi ■ d/dqi, velocities and Lagrangian transform locally as 



5qi 



dSqi 
"dT 



dSt 
'~dt 



SC = est + {dCldqi)5q^ + {dC/dqi) 



dSqi 



dSt 



dG ^ d{St) _ 
'dt' dt + ' ■ ^ * ~ 



dt " dt 

in terms of the total derivative of the generator of the transformation, the Noether charge 

G := C ■ 5t + pi ■ {5qi - iiSt) , 



(1) 



(2) 



and the Euler-Lagrange variational derivative Vi 
end points is 



dC/dqi — d(dC/dqi)/dt. The variation in action between fixed 



(55i2 = ^ dtSC^J^ dt 



dG ^ d{St) _ 
— -C-^+V,-{Sq^~q,dt) 
dt dt 



G(l)-G(2)- 



dt 



Sqi-Vi 



fdn dc\ 



(3) 



after integrating the term in d{St)/dt by parts. The action principle asserts that this variation vanishes for independent 
variations 6qi, St that vanish at the end points. It implies the Euler-Lagrange equations Vi = and dT-L/dt = —dC/dt, 
the rate of change of the Hamiltonian in non-holonomic systems. 
On-shell, where Vi = 0, 



<55i2 = ^ ~5Cdt = G(l) - G(2) 



dG 
'dt 



= 6C ■.= SC + C-{d6t/dt). 



(4) 
(5) 
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This extension of Noether's Theorem describes the evolution of any generator or Noether charge, in terms of the 
transformation of the Lagrangian it generates. It expresses the Euler-Lagrange equations of motion as the divergence 
of a Noether charge, which vanishes for a variational symmetry, but not for any other symmetry transformation 
(Section I). 

Our primary purpose is to contrast these different ways of reducing second-order differential equations to first-order, 
by comparing two familiar physical examples: 

Central Field Motion in a Static Potential which is completely integrable by virtue of energy and angular mo- 
mentum conservation, whether or not the system is scale invariant (Section II). 

Hydrostatic Gaseous Spheres in Adiabatic Equilibrium which are integrable only if they are scale invariant 
(polytropes) (Section III). 

Although generally not a conservation law, any symmetry of the equations of motion implies a useful dynamical 
or structural first-order equation Scaling symmetry is the most general simplification that one can make for 
any dynamical system. For the radial scaling transformations we will consider, Sr = r, the Lagrangian scales as 
some scalar density dC — —2uj£ and the action scales as SS = {1 — 2ui)S. The Noether charge generating the scale 
transformation evolves according to a non-conservation law ^ = (1 ^ 2Cj)C H, a first-order equation encapsulating 
all of the consequences of scaling symmetry. All the published properties 043 of index-n polytropes follow directly 
from this first-order equation. 

Our secondary purpose is to present our original variational formulation of spherical hydrostatics, our direct ap- 
plication of the scaling non-conservation law deriving from our extension of Noether's Theorem to non-variational 
scaling symmetry, our definition of a core radius, inside which all polytropes exhibit a common mass density structure 
(Sections III, IV). Section V completes the integration of the Lane-Emden equation by quadratures and obtains useful 
approximations to the Emden function 0„{^). 

An appendix reviews the thermodynamic properties of the physically important polytropes of index n = 3 @, |3i IE] ■ 
What is original here is the explanation of the the differences between relativistic degenerate white dwarf stars and 
ideal gas stars on the zero-age main sequence (ZAMS), following from their different entropy structures. Our original 
approximations to 03 (^) should prove useful in such stars. 

II. VARIATIONAL SYMMETRIES IMPLY CONSERVATION LAWS 

Time translation and spatial rotations are variational symmetries of the action integral, 

C{r,r,e)dt , C:=T{r,r,6)-V{r)^j{r^+r^0^)-V{r) , (6) 

so that the energy and angular momentum 

E{r,r)^{m/2)r^ + V{r) = {m/2){r^ + {l/mrf) + V{r) , I = mr^6, (7) 

are conserved. Because of these two first integrals, conservative central field motion is completely integrable by 
quadrature 



0{r)=da+ dr/{^2mr'^[E - V{r)]/P - r^} (8) 



t{r)^to+ dr/y/2r^[E - V{r)] - [l/mf, (9) 



where ''o are initial values at time t^. 

These two first integrals imply the first-order differential equations 



l/mr\ f^\l-[E~V{r)]~{—Y. (10) 



and the orbit equation 

dr 



^„ V'2mr4[£;-y(r)]//2-r2. (11) 
du 
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TABLE I: Period- Amplitude Relations and Virial Theorems for Inverse Power-Law Potentials V ~ 1/r' 



Tl 




Pprind-fnnnlitiiHp rplfltinn f ^l^^/*^ 


Virial theorem 


-2 


isotropic harmonic oscillator 


period independent of amplitude 


{K) = 


(V) 


-1 


uniform gravitational field 


falling from rest, e.g., z — gt^ /2 




{V)/2 





free particles 


constant velocity r ^ t 


{K) = 





1 


Newtonian potential 


Kepler's Third Law ~ 




~{V}/2 


2 


inverse-cube force 






-{V) 



What additional consequences follow if, V{r) ~ l/r", so that the equations of motion are also invariant under the 
infinitesimal scale transformation, 

5t={l + n/2)t , Sy^V; 6r^{~n/2)r, (12) 

which is not a variational symmetry of the action? Instead of another conservation law, scale invariance implies 
i/7.(i+"/2) — constant and the period-amplitude relations in Table L Because the kinetic and potential energies 
transform infinitesimally as 

5K = -nK , SV = -nV, (13) 

the time derivative of the virial A ^ pi • of a many-body system obeys 

A^2K + nV . (14) 

In a bounded system its time average < v4 = >, so that the time averages < K >, < V > obey the generalized 
virial theorems in the last column of Table I ^] . 



III. SCALING SYMMETRY MAKES HYDROSTATIC STELLAR STRUCTURE INTEGRABLE 

A. Variational Principle for Spherical Hydrostatics 

A non-rotating gaseous sphere in hydrostatic equilibrium obeys the equations of hydrostatic equilibrium and mass 
continuity 

— dP/ pdr — Gm/r'^ , dm/ dr = Ai:r^ p , (15) 

where the pressure, mass density, and included mass P{r), p{r), m(r) depend on radius r. As dependent variables, 
we prefer to use the gravitational potential V{r) = Gm/r'^dr and specific enthalpy (ejection energy, thermostatic 

potential) H{r) — J^''^^ dP/p, so that (fTS]) and its integrated form 

- dH/dr = dV/dr, V{r) + H{r) = , (16) 

R 

express the conservation of the specific energy as sum of gravitational and internal energies, in a star of mass M 
and radius R. In terms of the enthalpy H{r), these two first-order equations ([T5|) are equivalent to the second-order 
equation of hydrostatic equilibrium p7|) 

which is Poisson's Law for the gravitational potential V{r) = —H{r) — 

Because p{r), P{r), H(r) are even functions of the radius r, at the origin to order r^, spherical symmetry dP/dr = 
and mass continuity requires, 

p{r)^p,{l-Ar'), m(r) = ^.(l-^^r2) = l!^.p2/5^(r)3/5 . (18) 

3 5 3 

Thus, the average mass density inside radius r is p(r) := — pV^ p{rY^^ ■ 
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The second-order equation of hydrostatic equihbrium (jl7|) follows from the variational principle 5W = minimizing 
the Gibbs free energy 

W:= drC{r,H,H') (19) 

^0 

■ The Lagrangian 

C{r,H,H') = 4TTr^[-H'^/8TTG + P{p)] , '-.^d/dr , (20) 

is the sum of the gravitational and internal specific energies per radial shell dr. The canonical momentum and 
Hamiltonian are 

m dC/dH' = -r'^H'/G , Tiir, H, m) = ~Gm^/2r^ - ATir'^P{H) (21) 

arc the included mass and energy per mass shell. The canonical equations are 

dn/dm = H' = -Gm/r'^ , dU/dH = -m' = -iirr'^p . (22) 

Spherical geometry makes the system nonautonomous, so that dH/dr ~ —dC/dr = —2C/r vanishes only asymptoti- 
cally, as the mass shells approach planarity. 

B. First-Order Equation Between Scale Invariants 

The equations of hydrostatic equilibrium (15) can always be rewritten 

d\ogu/d\ogr ~ 3 — u{r) — w{r) , d\ogw/d\ogr = u{r) — 1 + w{r)/n{r) — d\og[l + n{r)]/d\ogr , (23) 
in terms of the logarithmic derivatives 

u(r) := dlogm/dlogr, v{r) := —d\og{P/ p)/d\ogr, w{r) := n{r)v{r) — ~d\og p/d\ogr , (24) 
and an index n{r) 

n{r):=d\ogp/d\og{Plp) , 1 + J- dlogP/dlogp , (25) 

n[r) 

which depends on the local thermal structure. Our homology mass density invariant w{r) will make explicit the 
universal mass density structure of all stellar cores, which is not apparent in the conventional pressure invariant w„. 

Because stars never have uniform mass density {n(r) ^ 0), their action cannot be invariant under radial translation. 
A hydrostatic structure will still be completely integrable, if the structural equations (15) are invariant under the 
infinitesimal scaling transformation 

5r ^ r, Sp = -nConP, SH = -CjnH, 5H' = -(1 -I- w„)ff', := 2/(n - 1) , (26) 

generated by the Noether charge 

7:r/2 prrrr/ 

G^:=-n-r-m- {(L„H) ^ r%— + 4nP{H)) ■ r + ^^^]- (27) 

The Lagrangian (20) then transforms as a scalar density of weight — 2cZ;„ 

SC = -2Cjn£ , SSi2 = (1 - 2Cjn) ■ Si2 , (28) 

so that only for a)„ = 1/2, the n = 5 polytrope , is the action invariant and scaling a variational symmetry. 

Both structural equations (|23p will be autonomous, if and only if n is constant, so that, P{r) = Kp{rY^~ , with 
the same constant K (related to the entropy) at each radius. When this is so, 

du/ dlogr — uii — u — Wn) , dwn/ dlogr — Wn{u ~ 1 -\ ) (29) 

n 

dlogwn dlogu dlogm 

— — diogr = (30) 



(3 - M - w„ ) 
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FIG. 1: Polytrope density gradient steepens as the boundary is approached (it — >■ 0). All solutions are tangent to the same 
density structure ■w„{z) ws = (5/3) (3 — u) at the center {u — 3), but differ for u < 2, outside the core radii marked by 
red dots. Approaching the outer boundary (it — )■ 0), the density pn{r) falls rapidly, so that its gradient Wn n[oaj"~^/u]^''"' 
diverges, for n < 5. 



0. 

In this section, we consider only the first equality (30), the first-order equation 

dWn Wn{u ~ 1 + Wn/n) 



(31) 



du m(3 — U — Wn) ' 

encapsulating all the effects of scale invariance We will consider only polytropes with finite central density pc, 
so that the regularity condition (fT8|) requires that all Wn{u) be tangent to |(3 — u). Such Emden polytropes are the 
regular solutions Wn{u) of the first-order equation, for which w„(u) — |(3 — u) for u — )■ 3. In place of u, we now 
introduce a new homology invariant z := 3 — u = —dlogpn/dlogr, where p„ := 3m{r)/4Trr^ is the average mass 
density inside radius r. In term of z, Wn, the characteristic differential equations (|30p are 

dz dlogWn dlogm 

= dlogr = — . (32) 



{3-z)iwn~z) 2^z+^ ^ 3 

Incorporating the boundary condition, the first of equations (32) takes the form of a Volterra integral equation Q 

Wn{z) - dz Wn^^^T^^!^^ - (5/J„)[l - (1 - Z/Sy-] WnPiciz) , Jn (971 - 10)/(7 - 71). (33) 

Jo (3-z)(w„-z) 

On the right side, the Picard approximation is defined by inserting the core values u;„(z) « (5/3)2; inside the integral. 
For ri = 0, 5, this Picard approximation is everywhere exact. For intermediate polytropic indices < ri < 5, the 
Picard approximation breaks down approaching the boundary, where Wn diverges as w„ — > n[ow"~^/M]^/", and is 
poorest for ?i « 3. Figure 1 shows the exact Wn{u) for n=0, 1, 2, 3, 4. 5. 

The second-order equation of hydrostatic equilibrium (|17|) . takes the dimensionless form of the Lane- Emden equation 

A.{e'^),eK-o . (34) 

in terms of the dimensional constant, dimensional radius and dimensional enthalpy 

- %^^Py""' > ^ := r/a , H = H,0n , where H, = {n + 1)— = (71+ l)ifpj/" . (35) 
47rG Pc 

The dimensionless enthalpy is On{^) and the dimensional radius, central density, included mass, mass density, average 
included mass density, specific gravitational force are 

771 17' 1 

r p„ mir)^ATTpca'-{-een), P»(r) = Pc-C(0, Pn{r) := = p,- (-3^/0, gir) ^ iirpca'i-e'J (36) 
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TABLE II: Scaling Exponents, Core Parameters, Surface Parameters, and Mass-Radius Relations for Polytropes of Increasing 
Mass Concentration. Columns 3-5 are well-known but columns 6-7 present a new measure of core concentration. 



n 


w™ 




Pcn{R)/PniR) 






T^ncore j 


Radius-Mass Relation ~ "/o^^n 





-2 


2.449 


1 


0.333 


1 


1 


7? ~ M^^^\ mass uniformly distributed 


1 


zilCX) 


3.142 


3.290 




0.66 


0.60 


R independent of 


1.5 


4 


3.654 


5.991 


132.4 


0.55 


0.51 


R ~ M~^/3 


2 


2 


4.353 


11.403 


10.50 


0.41 


0.41 




3 


1 


6.897 


54.183 


2.018 


0.24 


0.31 


M independent of R 


4 


2/3 


14.972 


622.408 


0.729 


0.13 


0.24 




4.5 


4/7 


31.836 


6189.47 


0.394 


0.08 


0.22 




5 


1/2 


oo 


oo 








0.19 


R — oo for any Af;mass infinitely concentrated; 



where prime designates the derivative ' :— d/d^. The scale invariants are 



The Noether charge 



G„(0 



711+1 



H2 n+l/ S 



evolves radially according to 

dGn _ _Hl d_ 
~df ~ ~~G ' d4 



1' 



qn+l 



1 



(37) 



(38) 



(39) 



This non-conservation law expresses the radial evolution of energy density per mass shell, from entirely internal 
(6',"+^/(n -I- 1)) at the center, to entirely gravitational (0^/2) at the stellar surface. 
For n = 5, scaling is a variational symmetry and (39) reduces to a conservation law for the Noether charge 



Gs = ^ • em^-f + f ) + = § • i^vlf^' ■ [-^5 - u/3 + 1] = constant 



(40) 



For the Emden solution, V5 is finite at the stellar boundary u = 0, the constant vanishes, so that W5{u) = 5w5 = |(3— w) 
everywhere. 

For n < 5, v„ diverges at the stellar radius ^1 , but lo„ — )■ oujn, a finite constant characterizing each Emden function. 
At the boundary u — 0, our density invariant w„(u) diverges as n[ouj 

n-i/y]i/n^ and 



(41) 



Table II lists these constants ^1, ow„, along with the global mass density ratios Pcn{R) / Pn{R) = { ^^^ )i and the 
ensuing dimensional radius-mass relation A/"^^ — [G/{n + 1)K]'^ ■ (47r/oWn"^^)-R"^'^ ■ Figure 1 shows Wn{u) for 
ri = 0, 1, 2, 3, 4, 5. Besides the well-known 0-13 third, fourth and fifth columns referring to the surface, we have added 
the sixth and seventh columns referring to the core radius, to be discussed in Section V. We have calculated all of 
Table II directly from the first-order equation (jSTJ, encapsulating all the effects of scale invariance j2l 



IV. EMDEN SOLUTIONS AND THEIR APPROXIMATIONS 



After obtaining Wn{z) :— —dlogpn/dlogr, either numerically or by Picard approximation, another integration gives 

I 

= »- {- r K(')-:p-.) } - f"^' 

0. - ^ »P {- 1' .„|„,„,'1)!;',(L.) } - <i - '/'f'" - (43) 
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TABLE III: Taylor Series and Picard Approximations 6nPic((.) to Emden Functions On{£,) 



n 


Emden Function 9n(£,) and Taylor Series 


Af„ := 5/(3n - 5) 


Picard Approximation 6'„pic(0 (1 + CVSA^n) 







-1 


1-CV6 


1 


sin e/C = 1 - ?V6 + CVl20 - CV5040 + ■ • ■ 


-5/2 


(1 - ^Vis)^''^ = 1 - CV6 + ?Vi2o - eViosoo + ■ • ■ 


n 


1 - CV6 + nC''/120 - n(8n - 5)/15120C'' + ■ ■ • 


5/(3n - 5) 


(1 + CVeA^n)"^" = 1 - + <Vl20 - n(6n - 5)CVl0800 + 


5 


(i + CVs)-'/' 


1/2 


(i + CVs)-''" 



m(z)/Af = C-f/' ■ exp{ /'d^^Tl 1 " f}} " (^4) 

3 ^Jz V[wn{z)-z\ 2z)) 3 

r(z)/i? = e/an-(f)^/^-exp{ fdzj \ 1 }U (|f)^ . (45) 

3 L (3 — z)[w„(z) — zj 2zJ J 3 — z 

All the scale dependance now appears in the integration constants M and R{M), which depends on M, except for 
n — 3. 

The Picard approximations 

0nPic(O = (l+e'WO^"^" , 7V„ :=5/(3n-5) (46) 

to the Emden functions are defined by inserting the core values Wniz) ss (5/3)z inside the integral and tabulated in 
the last column of Table IIL For polytropic indices n = 0, 5, this closed form is exact. For intermediate polytropic 
indices < n < 5, the Picard approximation remains a good approximation through order but breaks down 
approaching the outer boundary. Unfortunately, the Picard approximation is poorest near n = 3, the astrophysically 
most important polytrope. 

This n — 3 polytrope, which is realized in relativistic degenerate white dwarfs and in the Eddington standard 
model for luminous zero-age main sequence (ZAMS) stars, is distinguished by a unique M — R relation: the mass 
M = (-\/47r/oW3)(G/4i^)'^/^ is independent of radius R, depending only on the constant K := P j p^l"^ . In these stars, 
the gravitational and internal energies cancel, making the total energy = 17 + J7 = 0, leaving them in neutral 
mechanical equilibrium at any radius. 

Figure 2 compares three approximations to this most important Emden function, shown in yellow, whose Taylor 
series expansion is 

^3(0 = 1 - ^V6 + CV40 - (19/5040)^^ + (619/1088640)^^ - (2743/39916800)^^" + • • • . (47) 
Tenth-order polynomial approximation to this Taylor series expansion 

1 - 0.1666667^^ -h 0.025^^ - 0.0037698^^ ^ 0.0005686^* - 0.00006872^^° , (48) 
shown in red, diverges badly for ^ > 2.5 « 1.7 ^score. 
Picard approximation 

6'3Pic(e) - (1 + 2eVl5)"^/* = 1 - + ^V40 - 13^73600 • • • (49) 
1 - 0.1666667^^ 0.025^"* - 0.003611^^ -I- • • • , (50) 

shown in dashed green, converges and remains a good approximation over the bulk of the star, with < 10% 
error out to ^ ss 3.9, more than twice the core radius and more than half-way out to the stellar boundary at 
^13 = 6.897. This approximation suffices in white dwarf and ZAMS stars, except for their outer envelopes, 
which contain little mass and are never polytropic. Because it satisfies the central boundary condition, but not 
the outer boundary condition, the Picard approximation underestimates and overestimates 6'(^) outside 

C~3.9. 

Pade rational approximation Q: 

03Pad = ^1 ~ ^ll^Ti, l^^l?finn? = l-0.166667^2^0.025e^-0.00376984e^-|-0.0005686C^-0.0000857618ei°+- • • , 
1 + 17^^/108 + 4^/1008 

(51) 

shown in dashed heavy black, is a much better and simpler approximation. In fact, this Pade approximation is 
almost exact out to — 6.921, very close to the outer boundary ^13 — 6.897. 



These simple analytic approximations to ^3(^)1 simplify structural modeling of massive white dwarfs and ZAMS stars. 
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FIG. 2: The exact Emden function 63 (^) (solid yellow) and its polynomial (red), Picard (green dashed) and Pade (solid black 
dashed) approximations. Even in this worst case, the Picard approximation is holds out to twice the core radius at 2^3coi e = 3.3, 
before breaking down near the boundary. The Pade approximation is indistinguishable from the exact solution, vanishing at 
^1 — 6.921, very close to the exact zero ^13 — 6.897. 



V. A NEW MEASURE FOR CONCENTRATION INCREASING WITH POLYTROPIC INDEX 

Emden functions are the normalized solutions of the Lane-Emden equation p4p for which the mass density is finite 
at the origin, so that 0„(O) — 1, 0'^{0) = 0. Each Emden function of index n is characterized by its first zero On{^in), 
at dimensionless boundary radius As a new measure of core concentration, we also define the core radius ^core 
implicitly by M(^core) '■— 2, where gravitational and pressure gradient forces are maximal and w„ ~ 2, Pncore/Pnc ~ 0.4 
for all polytropes n > 1. Inside the core, the specific internal energy density dominates over the gravitational potential, 
so that for n > 1, 

u;„(u) « wsN = , 6'„(e)«l-eV6 , for u > 2, C < Ccore , (52) 

consistent with the universal density structure (18) all stars enjoy near their center. Inside the core, the enthalpy 
H{r) ~ On{£,) decreases as exp (— C^/6) for all polytropes. Outside the core, the gravitational potential V(r) dominates 
as it increases towards —GM/R at the stellar surface. 

The core concentration is illustrated in Figures 3 and 4, which show the local density pn/ Pnc = ^" as function of 

included mass fraction = (^2g")^ and of fractional radius -j^ — respectively. On each curve n in Figures 1, 3, 4, 
the core radii is marked by a red dot. 

The sixth and seventh columns in Table II list dimensionless values for the fractional core radius rncore/ R — £,ncore/£,i 
and fractional included mass mncore/M . 

For n 0, the mass is uniformly distributed, and the entire star is core. 

As < rt < 5 increases, the fractional core radius shrinks rcore/R 0, the mass concentrates towards the center, 

nincore/M ^ 0.19: 

For 1 < n < 3, the radius R decreases with mass M. Nonrelativistic degenerate stars have n = 3/2. 

For n=3, the radius R is independent of mass M. This astrophysically important case will be discussed in 
Section V and the Appendix. 

For 71 > 3, the radius R increases with mass M. As n — > 5, the stellar radius increases <^i„ -> 3(7i + l)/(5 — n), 
the core radius shrinks £,core ■\/10/3n, the fractional core radius rcore/R = ^core/£,in ~^ 0.045(5 — n), 
rUncore/M 0.19, and o^n v^3/^i„ 0. 
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m/M 



FIG. 3: Normalized density profiles as a function of fractional included mass m/M, for polytropes of mass concentration 
increasing with n. The red dots mark the core radii, at which the densities stay near p{rcore) / pc ~ 0.4, for all n > 1. For 
uniformly distributed mass (n = 0), the polytrope is all core. As the mass concentration increases (n — )■ 5), the core shrinks to 
20% of the mass. 

For n = 5, the stellar radius R — oo for any mass M, while the included mass m(r) is concentrated towards finite 
r. Scaling becomes a variational symmetry, so that the Noether charge G5 in (36) is constant. For the Emden 
solution, this constant vanishes 

G,-m^f + ^-^) + lerj',]^^e,e',{v,-u/3-i)^o , (53) 

SO that t;5 = l — u/3, 0'^ = — Integrating then yields 

05io = a+e/3)-'^' : P5 = P5c{i+em-'/' , (54) 

after normalizing to ^5(0) = 1. 
For n > 5 the central density diverges, so that the total mass M would be infinite. 

VI. CONCLUSIONS 

We have distinguished two ways in which symmetry reduces a second-order differential equation to a first-order 
equation: 

Variational Symmetry (epitomized by central field dynamics): The symmetries of the action lead to conser- 
vation laws, first integrals of the original equations of motion 
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0.0 0.2 0.4 0.6 0.8 1.0 



r/R 

FIG. 4: Normalized density profiles as function of fractional radius rjR. The density is uniform for n = 0, but is maximally 
concentrated at finite radius for the n = 5 polytrope, which is unbounded (R — oo). The density at the core radius stays about 
p{rcore)/pc ~ 0.4, for any n > 1. 

Non-variational Symmetry of the Equations of Motion (epitomized by scale invariant hydrostatics): 

Yields a first-order equation between scale invariants which, although not a first integral, still leads to 
integration by quadratures. 

In the latter case, we obtained all the familiar properties of polytropes, directly from the first-order equation between 
invariants. We observed that, like all stars, polytropes of index n share a common core density profile and defined a 
core radius outside of which their envelopes differ. The Emden functions On{^), solutions of the Lane-Emden equation 
that are regular at the origin, were finally obtained, along with useful approximations. 

The Appendix reviews the astrophysically most important n = 3 polytrope, describing relativistic white dwarf stars 
and zero age main sequence stars. While reviewing these well-known applications 15|, we stress how these same 
mechanical structures differ thermodynamically and the usefulness of our original (Section V) approximations to these 
Emden functions. 



APPENDIX: IMPORTANT ASTROPHYSICAL APPLICATIONS OF n = 3 POLYTROPES 

The n = 3 polytrope, which is realized in white dwarfs of nearly maximum mass and in the Eddington standard 
model for luminous zero-age main sequence (ZAMS) stars that have just started to burn hydrogen, is distinguished 
by a unique M — R relation: the mass Af = 47r qlo^ (K j ■kG)^!'^ is independent of radius i?, depending on the constant 
K := Pj p^l'^ . In these stars, the gravitational and internal energies cancel, making the total energy = + C/ = 0, 
leaving them in dynamical equilibrium at any radius 0] • 

1. Relativistic Degenerate Stars: K Fixed by Fundamental Constants 

White dwarfs of nearly maximal mass are supported by the degeneracy pressure of relativistic electrons, with number 
density rie = p/^efnu, where mn is the atomic mass unit and the number of electrons per atom /ie = Z/A — 2, 
because these white dwarfs are composed of pure He or C^^/O^^ mixtures. Consequently, Kwd = ^[f-]^^^™HA'e"'*^'^ 
depends only on fundamental constants. This universal value of Kwd leads to the limiting Chandrasekhar mass 
= ^Mjf^l = 5.824Mo//.2 = l.456Mo • (|)2 |,[3. 
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2. Zero-Age Main Sequence Stars: K{M) Depends on Specific Radiation Entropy 

In an ideal gas supported by both gas pressure Pgas = TZpT/fj, :— /3P and radiation pressure Prad = :— 
(1 — I3)P, the radiation/gas pressure ratio is 

Pgas P p Sn 

The specific radiation and ideal monatomic gas entropies are 



4aT3 n , T{rf' ^ 



^rad = ^, 5gas(r) =(-)-log(^^) , (A.56) 



so that the gas entropy gradient 



depends on the difference between the adiabatic gradient V ad = 2/5 and the star's actual thermal gradient V := 
dlogT/dlogP, which depends on the radiation transport. 

Bound in a polytrope of order n , the ideal gas thermal gradient and gas entropy gradient are 

For n > 3/2, the thermal gradient is subadiabatic, the star's entropy increases outwards, so that the star is stable 
against convection. 

Zero-age main sequence stars (ZAMS), with mass OAMq < M < 15OM0, have nearly constant radiation entropy 

S'rad (Af), because radiation transport leaves the luminosity generated by interior nuclear burning everywhere propor- 
tional to the local transparency (inverse opacity) k^^. Assuming constant Siadl-^^)? we have Eddington's standard 
model, an n = 3 polytrope with Si-^^iM) = 4(7?.//i) • (1 — /?)//? and 

K{M) = P/p^'^ = {[3(1 - f5)/a]{^l/^iPfY'^ (A.59) 

depends only on (3{M), which is itself determined by Eddington's quartic equation 043] 



1-/3 fMp^Y nr 3v;iO^/_?i^x3/2 



The luminosity L = LEdd[^ ~ depends on the Eddington luminosity LEdd '■— 4:t:cGM/ Kp through the 

photospheric opacity Kp. From Eddington's quartic formula, the stellar luminosity 

L = L EddiO -003) p^lS^M/Mef. (A.61) 

This mass-luminosity relation is confirmed Q in ZAMS stars: On the lower-mass ZAMS, /3 w 1, L ^ M^; on the 
upper-mass ZAMS, /3 « (^) M. 
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